Inducing spin-correlations and entanglement in a double quantum dot through non-equilibrium 
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For a double quantum dot system in a parallel geometry, we demonstrate that by combining the effects of 
a flux and driving an electrical current through the structure, the spin correlations between electrons localized 
in the dots can be controlled at will. In particular, a current can induce spin correlations even if the spins are 
uncorrelated in the initial equilibrium state. Therefore, we are able to engineer an entangled state in this double- 
dot structure. We take many-body correlations fully into account by simulating the real-time dynamics using the 
time-dependent density matrix renormalization group method. Using a canonical transformation, we provide an 
intuitive explanation for our results, related to Ruderman-Kittel-Kasuya-Yoshida physics driven by the bias. 
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Introduction: Considerable progress in nanotechnology in 
the last decades has made possible the fabrication of new ar- 
tificial structures such as quantum dots (QDs), quantum 
rings, or molecular conductors. The physics of quantum dots 
in a parallel geometry is intriguing, since it allows one to 
study interference effects between electrons traveling through 
different paths, most notably realized in the Aharanov-Bohm 
effect. Such structures have been studied in several experi- 
ments fU,^. Besides the interest in practical applications in 
nanoelectronics or in fundamental many -body physics such as 
the Kondo effect [3], double quantum dots (DQD) also play a 
vital role in the context of quantum information processing 
Generating, controlling, and detecting entangled states 
in condensed matter systems is one of the challenges for fu- 
ture quantum computation applications 10]. Various proposals 
for entangling spatially separated electrons have been put for- 
ward, such as, for instance, by splitting Cooper pairs Q8l |9[ or 
by manipulating spins in quantum dots jlfJl , In a DQD, an 
entangled state can be realized by putting the electrons into a 
singlet state flE, 11]. Means of detecting entangled states of 
electrons were discussed in, e.g., Refs. 



In this work, we demonstrate that an entangled state 
between electrons localized in a DQD embedded in an 
Aharonov-Bohm interferometer can be induced and con- 
trolled by sending an electrical current through the struc- 
ture. In the presence of a flux, the initial state can even 
be fully uncorrelated yet the non-equilibrium dynamics re- 
sults in non-zero spin correlations in the steady state. The 
sign and the strength of such steady-state spin correlations 
depend on voltage, interactions, and the flux. The genera- 
tion of entanglement through non-equilibrium dynamics in 
quantum dots, with different set-ups, has been discussed in 
Refs. Ill 311 . An additional motivation for our work stems from 
the current interestjboth from theory lll4l U5I and from exper- 
iment (see, e.g., 11611 ). in the real-time and non-equilibrium 
dynamics of nano-structures with strong electronic correla- 
tions. An important aspect is that we can treat both interac- 
tions and non-equilibrium dynamics in a well-controlled man- 
ner using the time-dependent density matrix renormalization 
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FIG. 1 . (Color online) Current and spin correlations as a function 
of time, (a) Current J(t) for (f> = and it, V = 0.50. (b) Spin 
correlation Si 2 (t). All for N = 41,U = 0.5, V = 0.5. Inset in (b): 
Sketch of the DQD structure studied here. 



group (DMRG) method Jnll . 

The model: We study a system of two quantum dots in a 
parallel geometry represented by two Anderson impurities, as 
sketched in Fig.QJb). The Hamiltonian is 
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The system size is 2N + 2 where TV is the number of sites 
in the left or right lead. The two dots are at the center of 
the system labeled by j = 1, 2. The Hamiltonian consists of 
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four parts: First, the non-interacting leads Hi with a constant 
hopping matrix element to = 1 used as the unit of energy 
(% = 1, e = 1). Second, the terms -ffhyi an d -^hy2 gi ye ri se to 
the hybridization between the localized levels of the dots and 
the leads. If not stated otherwise, we consider fully symmetric 
tunnel couplings, i.e., t\ = t' 2 = t' (see the supplementary 
material fll 811 for a discussion of asymmetric couplings). We 
define the tunneling strength by Y = 2irt' 2 pi ea d s {E f) = 2t 12 , 
where pi ads(^_F) is the local density of states (LDOS) of the 
leads at the Fermi energy Ep. In the hopping matrix element 
between the second dot and the right lead we incorporate an 
arbitrary phase <f> that will be important in our study. Finally, 
there is the interacting region H lnt with the two quantum dots, 
which are both subject to the same Coulomb repulsion U and 
a gate potential V g = —U/2 such that both dots are kept at 
half filling. The operator c} alo (c a i a ) creates (annihilates) an 
electron at site I in the a = L,R lead with spin a while 

(d ja ) acts on dot j; n aia = c f al<j c a i a as usual. In Eq. ©, [i L 
and [in mimic the chemical potentials of the leads. 

The ground state and the linear conductance of DQDs 
Eq. (Q3 were extensively studied in Ref. Ill 911 . A closely re- 
lated DQD model with a finite flux <fi and with spin-polarized 
electrons was discussed in Ref. l20ll . 

The phase included in Eq. (|4]i may have a different meaning 
depending on the specific physical realization. The most obvi- 
ous one is to associate <fi with a magnetic flux that pierces the 
ring structure containing the two dots and the first site from 
each lead as shown in Fig. Q2b). As usual, one can use a 
gauge transformation such that the flux appears in only one 
of the four hopping matrix elements. 

Another situation that can be described by this Hamiltonian 
is a single quantum dot with two levels. Depending on the fill- 
ing of a QD, two degenerate levels may exist that can couple 
differently to the leads. As an example of this latter case let 
us consider a spherical QD weakly coupled to the leads. The 
angular part of the QD's eigenfunctions are spherical harmon- 
ics Yim- Let us suppose that the weak contact with the leads 
are along the z-axis, and discuss the coupling of p^-orbitals 
(/j = x, y, z). The p z orbitals have an opposite weight for 
opposite z-coordinates creating a phase difference of <f> = n 
in the hopping terms to the leads. The other two orbitals p x 
and p y lack this extra phase by symmetry. This example il- 
lustrates that different phases for different orbitals of a single 
QD are possible once they are coupled to the leads through 
the hybridization. Note that in the first example, <f) can take 
arbitrary values (j> 6 [0,7r], whereas in the second case, only 
4> = 0, 7r are possible. 

We use DMRG lfl7tl to obtain the steady state in the pres- 
ence of a finite bias voltage by time-evolving the wave- 
function 1^(4)) and then its properties such as the current and 
spin correlations as a function of time t. This method has been 
successfully used to study non-equilibrium transport through 
nano-structures with electronic correlations lfl4l 2ll 22 1. We 
measure the spin correlations by evaluating H23II 
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FIG. 2. (Color online) Illustration of the canonical transformation 
used to understand the non-equilibrium spin correlations, (a) (j> = 0. 
(b) 4> = it- The upper two panels (a.l) and (b.l) show the situa- 
tion for V = 0, whereas (a.2) and (b.2) are for V 7^ 0. The ap- 
plication of the bias leads to a ladder structure where the bias acts 
like a transverse hopping matrix element between the symmetric (s- 
channel) and antisymmetric states (a-channel) defined in Eq. {8]l. 



The current between two sites in the leads is defined as 



Jl,m(t)=it J2(*( t )\ C l 



T ci a \V(t)). (7) 



Si 2 {t) = (*(t)|Si ■ Sa|*(*)) 



(6) 



In the figures, we display the current J = (Jli,li + 
Jri,R2)/2 averaged over the first link in the left and right lead. 

Our simulations start from the system in equilibrium with a 
finite T 7^ and a charge per spin of (rij a ) = 0.5 on both dots. 
At time t = 0, we turn on a bias voltage V = — M_r that 
drives the system out of equilibrium. We work at large values 
of T = 0.25 such that the transient dynamics to reach the 
steady state is short i22tl , During the time-evolution, we use 
time steps of St ~ 0.1 and enforce a fixed discarded weight 
] 3 of 10~ 5 or less, keeping a maximum of 2000 DMRG 
states. All runs are performed at an overall half filling of dots 
and leads. 

Results: In Fig. Q] we elucidate the time-dependence of 
the current and spin correlations, comparing the behavior of 
<fi = to (j> = it. Similar to a single quantum dot 12211 . the 
current undergoes transient dynamics, and then takes a quasi- 
stationary value (i.e., a plateau in time), which we shall refer 
to as the steady-state regime. Since we work with leads of 
a finite length, ultimately, the steady-state current decays and 
changes its sign (22]. The comparison with other methods 
shows that DMRG correctly captures properties of the true 
steady state as system size increases IU4I1 . 

For the spin correlations shown in Fig. Q3b), we first ob- 
serve that in the initial state, Si 2 > for <j> — whereas the 
correlation vanishes for <f> = n. The application of the bias 
voltage does virtually not affect the value of Si 2 for = 0, 
which remains positive in our simulations. The more interest- 
ing behavior is realized for <f> = it. As a function of time, S12 
decreases and approaches a roughly constant value. The tran- 
sient time seems to be comparable to the one for the current 
and is of order 1 /Y. Moreover, the transients are suppressed 
by increasing the bias, similar to a single quantum dot H22I1 . 
This finite and large spin correlation between the spins local- 
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FIG. 3. (Color online) Current and spin correlations as a function 
of bias V and phase 0. (a) (J) /V vs. V for = and (j> = ir 
in units of twice the conductance quantum Go- (b) Steady-state spin 
correlations (S12} vs. bias V. The figure shows data for different 
system sizes N = 25, 35, 41 for — ir, with very little finite-size 
dependencies. In (a) and (b), U/T — 2. Inset in (b): Steady-state 
spin correlations vs. U/T for V = 0.5. (c) Spin correlations (S12) 
vs. for V = and in the steady state with V = 0.5. 



ized in the dots that emerges in the steady state and that is 
induced by driving a current through the structure is the main 
aspect of our work. It implies that non-equilibrium dynamics 
can be used to prepare a DQD in a correlated and thus entan- 
gled state. 

The qualitative behavior of the spin correlations can be un- 
derstood by using a canonical transformation of the states of 
the leads, which is given by (see, e.g., |2^, 251): 



C 7 /a = (CRIct ± CUc 



(8) 



where 7 = s, a are the symmetric and antisymmetric com- 
binations, respectively. The result of this transformation is 
sketched in Fig. |2] where the leads shown there now repre- 
sent the new states obtained from Eq. ((H). In the absence of 
a bias voltage, there is no direct coupling between these new 
states, as depicted in Figs.[2|a. 1) and (b.l). Most importantly, 
the dots are coupled to only the symmetric states for = 0, 
whereas for <p = tt, dot j = 1 is coupled to the symmetric 
states and dot j = 2 to the antisymmetric ones. For = 0, the 
Ruderman-Kittel-Kasuya-Yoshida (RKKY) interaction gives 



rise to a ferromagnetic correlation between the dots since each 
path that connects them involves an odd number of sites and 
since the leads are at half-filling JI9I1 . For = tt, the dots are 
part of two decoupled subsystems and therefore, Si 2 vanishes. 

Upon applying a bias, one effectively obtains a ladder ge- 
ometry where the voltage acts as a transverse coupling be- 
tween the symmetric and antisymmetric states of Eq. (J8j as 
shown in Figs. |2ja.2) and (b.2). For = 0, the coupling 
V only marginally affects the correlations. By contrast, for 
= 7T and V 7^ 0, the dots are now connected through 
paths with an even number of sites in the effective leads and 
therefore, in the ground state of such a geometry, one ex- 
pects a finite negative spin correlation. Our numerical results 
shown in Fig. Q2b) unveil that the same behavior occurs in 
non-equilibrium as well. While here we focus on fully sym- 
metric tunnel couplings, the main results can be recovered in 
the case of asymmetric couplings fljll . and therefore, fine- 
tuning of parameters is not necessary to observe a change of 
Si 2 induced by a bias V. 

After qualitatively explaining the emergence of finite spin- 
correlations in the current-carrying stationary state, we next 
study the dependence of the steady-state properties on the 
bias potential. We denote the steady-state values by (S12) 
and (J), obtained from averaging over time-dependent data 
in the steady-state regime (compare Ref. 122J]). Figure |3a) 
shows (J) /V versus V for phases = and tt. For = 0, 



(J) /V approaches a constant value at low bias 111911 (an ex- 
tended analysis of the current-voltage characteristics of our 
DQD structure will be presented elsewhere). For = tt, the 
linear conductance vanishes due to the Aharanov-Bohm ef- 
fect lEoll . A finite voltage causes a finite current to flow in 
both cases, but (J) /V for = is always larger than in the 
= 7T case. 

In Fig. [3jb), we show the steady-state spin correlations 
(S12) versus V. First, let us emphasize that data for the 
steady-state values obtained from systems of different lengths 
are included, showing that all our main results are quantita- 
tively robust against finite-size effects. For = 0, a constant 
value of (S12) > is found. A slight decrease appears for 
V > 1, which we trace back to the variation of the LDOS 
of the leads seen by the dots. Since in our simulations we 
work with tight-binding bands with a finite band-width and 
band curvature, this LDOS decreases with V. For the case of 
= 7T, the value of the steady-state correlations can be tuned 
by the bias voltage. As the figure clearly shows, we can get to 
(S12) ~ —0.25 for U/T = 2. We therefore realize a mixed 
state with singlet correlations dominating over triplet corre- 
lations. The steady-state values further depend on U /T. To 
elucidate this, we plot (S12) versus U in the inset of Fig. |3jb) 
for a fixed value of V = 0.5. As expected, the larger U, the 
more strongly charge fluctuations are suppressed, leading to 
larger steady-state spin correlations. 

So far we have investigated the dependence of correlations 
on V and in non-equilibrium, comparing the cases of = 
to = n. An additional degree of tunability can be added 
if the phase can take arbitrary values. To illustrate this as- 
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FIG. 4. (Color online) Effective J G ff for DQD system as a func- 
tion of V (main panel) and U (inset) calculated for (j> = tt. Same 
parameters as in Fig.[3}b)- 
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pect, we display the static spin correlations (i.e., those at 

V = 0) and the steady-state correlations (here for V = 0.5) 
in Fig.[3]c). As expected from the discussion of Fig. El (S12) 
is positive for small <fr at V = and then decreases to zero as 
4> = 7r is approached. This transition to the uncorrelated case 
of = 7r is continuous. At a finite voltage, it is possible to 
go from positive steady-state correlations to negative ones by 
changing <f>. For the parameters of Fig. EJc), the steady-state 
correlations change sign at <fi c w O.I871". This value depends 
both on U and V. To summarize, the steady-state correlations 
can be tuned both in sign and magnitude by changing V, </>, 
and U. 

Based on the qualitative picture developed so far, we con- 
clude that the steady-state correlations are a result of mixing 
the symmetric and antisymmetric states of lead electrons in 
non-equilibrium. At finite U, this may be viewed as an RKKY 
effect in non-equilibrium. To estimate the strength J c g of the 
effective coupling between the localized spins in the dots in 
the steady state, we proceed as follows. First, in the set-up 
shown in Fig. |2jb. 1), we connect the impurities by an ad- 
ditional term J e gSi ■ S2 and we then calculate the equilib- 
rium correlations S12 = Si2(J e ff) in this reference system. 
In the next step, we use the numerically determined function 
(S12) = (Si2)(V) and by equating it to Si2(J e ff), we obtain 
Jett = Jes(V), keeping all other parameters fixed in both the 
time-dependent simulation and in the reference system. The 
results are presented in Fig.|4]and its inset. 

Finally, we study the robustness of the steady-state spin cor- 
relations against quenching certain parameters of the Hamil- 
tonian Eq. ((T). We proceed as before, i.e., a finite bias volt- 
age V > is turned on at t = 0, and in addition we perform 
quenches of either the phase <$> or some of the tunnel couplings 
at a time t q > 0. 

The first example, presented in Fig.[5]a), consists of chang- 
ing the phase from to tt. We show results for t q = and 10 
and observe that in both cases, the spin correlations change 
from their initial constant and positive value to a negative 
value that only depends on parameters in the steady state (i.e., 

V and <p), but as this comparison shows, is virtually indepen- 
dent of the specific transient dynamics. 

In the second case, presented in Fig. |3b), we simply dis- 



FIG. 5. (Color online) Behavior of the spin correlations in various 
quantum quenches, (a) Quench of the phase cj> from cj> = to <f> = 7T 
at a time t q — 0, 10 (dashed and dot-dashed curves), (b) Quench of 
the hopping matrix elements t[ and t' 2 to zero at t q = 20. This de- 
couples the DQQ from the leads, (c) Quench of just t[ to zero at time 
t q = 20, effectively decoupling dot 1 while the current continues to 
flow through dot 2. In all cases, we start with the usual protocol, i.e., 
we quench the bias from zero to V = 0.5 at t — 0. For comparison, 
the data from Fig.QJb) without any additional quenches are included 
(lines with open/solid circles are for <f> = n and = 0, respectively). 
In all panels, the time t q at which the quench is performed is indi- 
cated by arrows. In this figure, U/T = 2, V = 0.5, and TV = 41. 



connect the quantum dots from the leads at time t q > after 
the steady state has been established. As expected, the spins 
remain in a correlated state after isolating them from the re ser- 
voirs. In order to extract pure singlets from this mixed entan- 
gled state J26l |27m one would need to be able to reach spin 
correlations (Si 2) < —1/4 J27, 28], while charge fluctuations 
render the absolute value of steady-state correlations smaller 
in our case. 

Third, after reaching the steady state, we isolate one of the 
dots while the current continues to flow through the other one . 
This results in a destruction of the spin correlations after a 
short transient time [see Fig. 0c)]. Therefore, control over 
the tunneling matrix elements allows one to put the system 
back into its original uncorrelated state. Both the generation 
of entanglement and the removal happen on short times scales, 
similar to the proposals discussed in lfl3ll . 

Summary: In this work, we demonstrated that spin correla- 
tions between spatially separated electrons localized in a par- 
allel DQD embedded in the rings of an Aharonov-Bohm in- 
terferometer can be induced and modified by driving a current 
through the structure. The steady-state correlations depend on 
voltage, the flux, and Coulomb interactions. Control over the 
individual tunneling couplings would allow one to isolate the 
entangled spins from the environment or to remove the entan- 
glement again. The mechanism behind this time-dependent 
formation of correlations can be thought of as an RKKY effect 
in non-equilibrium. Our results may be relevant for applica- 
tions of DQD structures in quantum information processing. 
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Supplementary information for Inducing 
spin-correlations and entanglement in a double 
quantum dot through non-equilibrium transport 

ASYMMETRIC TUNNEL COUPLINGS 

In this supplementary material we show that the key results 
of our work described in the main text for a DQD with sym- 
metric tunnel couplings carries over to the asymmetric case. 
We start from the most general form of the Hamiltonian from 
Eqs. (3) and (4) of the main text. The hybridization between 
the localized levels of the dots and the leads can be written as 

-Hhyi = -t' L id[ a CLi<T - t' R1 d\ a c R i a + h.c, (9) 
H hy 2 = -t' L2 dla c Li<j - t' R2 e 1 ^ d\ a CR\a + h.c. , (10) 

where tunnel matrix elements t a j, are considered that can in 
general be different for each dot (j = 1 , 2) and lead (a = 
L, R). The tunneling strength is then defined by Tj = t'£j + 
t'jij, j = 1,2. For each dot we introduce an angle 6j via 
Oj = BXctB,n{t' L j/t' R j), and therefore, the terms from Eqs. (|9} 
and ( [Tol l can be reexpressed as 

Hhyj = d] a (sinflj c Lla + cos0j c R \ a + h.c.) . 

(11) 

We next generalize the canonical transformation of the lead 
operators that was discussed in the main text in the following 
way: 

c^ia = (sin 6»i c Lir7 + cos6»i c Ria ) , (12) 
Cvio- = (cos Ox c Li(T - sin 9i c Ria ) . (13) 

In the special case of 0i = 7r/4, this reduces to Eq. (8) from 
the main text with /i — s and v = a. Note that the oper- 
ator c^ia defined in Eq. (TTZb appears directly in -ffhyi from 
Eq. (fTTT i. Therefore, the QD j = 1 is just connected to the fx 
channel. The resulting configuration, i.e., QD 1 coupled to the 
/j-channel and QD j = 2 coupled to either the /i-, v-, or both 
channels is depicted in Fig. [6] 

The Hamiltonian of the leads H\ with chemical potential 
IXr = —/JtL = V/2 is, under this transformation, given by 

H l = H c + H v , (14) 

N-l 

H - = -*o E E ( c r- c 7'+i CT + h.c) (15) 

N 
i—l\a 

+^71^0-,] (16) 

with 

t p = -Vsindx cos 0i, (17) 
= -V v = ^ (cos 2 6x - sin 2 6^ . (18) 



-V/2 V/2 




FIG. 6. Sketch of the representation of the DQD structure and its 
coupling (a) the original leads and (b) the effective leads defined in 
Eqs. (12} and ( 1 1 3b - We refer to the effective leads as the fi- and v- 
channel, respectiveley. 

Similar to the case of fully symmetric couplings shown in 
Fig. 2 of the main text, the bias potential connects the chan- 
nels fj, and v through an effective coupling t p . This coupling 
takes its maximum value in the case of symmetric couplings 
forQD j = 1, i.e., t' L1 = t' m and therefore, 6\ = ir/4. 

Next, the second QD is connected to the new channels fi 
and v through: 

a 

with 

= (cos 6l 2 cos 0i + sin 2 sin0i e^) , (20) 

t'»2 = V^2 (-cos02 sin0x + sin0 2 cos0i c^) .(21) 

It is straightforward to verify that |^ 2 | 2 + \t' v2 \ 2 =T 2 . 

Figure [7] shows, for <f> = 0, tt/2 and it, the square modu- 
lus of t'^ and t' v2 as a function of 0i and 02 for U = 0.5, 
Ti = T 2 = 0.25, N = 35. For <j> = and 4> = tt we observe 
diagonal stripes where or t' v2 vanish resulting in the situ- 
ation depicted in Figs. 2 (a) and (b) from the main text. For 
(f> = it is obvious that, if 02 = 0i, the canonical transforma- 
tion will decouple both QDs from the z^-channel at the same 
time as both dots have the same asymmetry. 

In order to reproduce the features of Fig. 2 of the main text 
for the case of asymmetric couplings, we require that either 
t'fj.2 or t'v2 b e close t° zero. To illustrate this we have selected 
eight parameter sets for the tunneling matrix elements and the 
flux (j) that reproduce such situations, which we list in Table H] 
Each set of parameters contained in this table is labeled by a 
letter A,B, . . . , G. 

Set A is the completely symmetric situation with = that 
we discussed in the main text. In this case, t' 2 / ^/T^ = 1 and 
i' u2 = 0. Therefore, QD j = 1 and j = 2 are connected to the 
same effective lead and the spin correlations are positive in the 
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0.25 0.5 0.75 1 0.25 0.5 0.75 1 0.25 0.5 0.75 1 



FIG. 7. Square modulus of the couplings between QD j '■ — 2 and the effective channels [i and v defined in Eq. ( 112b and J 1 3 b - All curves 
are for T2 = 0.25, U = 0.5, V a = — U/2 and AT = 35. To reproduce the features of Figs. 2(a) and (b) of the main text, we need to require 
\t'v2\ 2 ~ or jt^l 2 — 0, respectively. Note that t'^ 2 and t' u2 do not depend on Y\. 



Set 


0l/7T 


02 A 


t'n 


t'rl 


i' 

l l2 


t'r2 




l^ 2 | 2 /r 2 


K 2 \ 2 /r 2 


t P /v 


Vy.IV 


(Sl2> 


A 


0.25 


0.25 


0.353 


0.353 


0.353 


0.353 


0.0 


1.000 


0.000 


0.500 


0.000 


> 


B 


0.09 


0.89 


0.139 


0.480 


0.169 


-0.470 


1.0 


0.996 


0.004 


0.267 


0.422 


> 


C 


0.48 


0.47 


0.031 


0.499 


0.497 


-0.047 


0.2 


0.997 


0.003 


0.067 


-0.496 


> 


D 


0.25 


0.25 


0.353 


0.353 


0.353 


0.353 


1.0 


0.000 


1.000 


0.500 


0.000 


< 


E 


0.14 


0.64 


0.213 


0.452 


-0.213 


0.452 


0.0 


0.000 


1.000 


0.385 


0.318 


< 


F 


0.07 


0.43 


0.110 


0.488 


0.488 


0.110 


0.5 


0.090 


0.910 


0.212 


0.452 


< 


G 


0.12 


0.63 


0.184 


0.465 


0.459 


-0.198 


0.2 


0.049 


0.951 


0.342 


0.364 


< 



TABLE I. Relationship between the bare parameters of the system [see Eqs. {9} and d 10H and the tunnel matrix elements of the effective 
channels \i and v [compare Eqs. Jilt and t!9H . In this table we show a few examples for which the discussion of the main article applies even 
for an asymmetric coupling. We choose Ti = Y 2 — 0.25 for simplicity. In the last column we list the expected sign of the spin correlations. 



steady state. A similar behavior is realized for sets B and C. 
Set B corresponds to a situation in which (p = tt while set for 
C, we choose <f> = Q.2n. Note that in set B the phase changes 
the sign of t' R2 and the situation is thus similar to set A. In set 

C, all four couplings, t' L1 , t' R1 , t' L2 , and t' R2 , are different from 
each other and for the selected value of the phase <j> = 0.2tt, 
no signs are affected. 

The examples labeled D to G are cases in which t' 2 is 
small, and therefore, where, according to the picture put for- 
ward in the discussion of Fig. 2(b) of the main tex, we expect 
negative spin correlations in the steady state. Set D is the com- 
pletely symmetric case with 6\ = 62 = 7r/4 and <fi = tt. This 
case was analyzed in the main article and is shown here for 
comparison. For set E, although = 0, there is a difference 
of tt/2 between 6\ and 82, resulting in a behavior similar to set 

D, despite the left-right asymmetry present in set E. Note that 
t' L2 is negative as a consequence, and therefore, the steady- 



state spin correlations are also negative, in accordance with 
the discussion of the main article. The set F realizes a case 
with <fi = tt/2 and a small t' u2 ~ 0.09. In this example, left- 
right symmetry is broken but there are still just two different 
values for the four tunnel matrix elements t' a ^ . Finally, in set 
G, <f> = 0.2tt and all four tunnel matrix elements are different. 
The very small value of t' u2 implies negative steady-state spin 
correlation as before. 

To verify our predictions for the sign of steady-state spin 
correlations that are based on the numbers listed Table H] we 
have calculated the spin correlation as a function of time. The 
results are shown in Fig. [8] Sets A, B and C, as they all have a 
small value of t' 2 , lead to a positive spin correlation positive 
in the initial state (t = 0) whose sign and absolute value are 
barely affected by the bias at all. Sets D, E, F and G all have a 
small value for t' u2 . As a consequence, the dots are predomi- 
nanntly connected to one of the two channels fi or v [compare 



8 



o.i - 




a* ax ax -ax ax axax ax a x .•i#,&:t#.&:dl.tx.'i»:&:&:ift,-df:tP:i!*:il*; i 



U = 0.5 . . A D 

u/r = u/r>2 t'Z*c -~-~f 



""1 —1 ~ x xC F 

V = 0.5 --- a 




o 



o 



5 



10 



15 
time 



20 



25 



30 



FIG. 8. Spin correlations as a function of time for the parameter 
sets presented in Table [I] Note that for sets with small values of t' V 2, 
such as A, B and C, the spin correlations do practically not change 
when the bias is applied. By contrast, whenever t^ 2 is small, as is 
the case in sets D, E, F and G, the spin correlations in the steady state 
are negative. Sets A and D are the cases discussed in the main text, 
compare Fig. 1(b). 

Fig. 2(b.l) of the main text] and spin correlations vanish in 
equilibrium, i.e., at for V = 0. Upon applying a bias, negative 
spin correlations emerge, as expected. Note the small value 
of the steady-state spin correlation in set F, consistent with the 
observation that in this case, the value of the coupling between 
the effective leads t p is the smallest. 

As a conclusion, even for asymmetric tunnel couplings be- 
tween the DQD and the leads the application of a voltage can 
induce and change spin correlations, which is the main result 
of our work. The main requirement to obtain a large effect of 
the bias V on the non-equilibrium spin correlations is a small 
value of either t^- 



